Inconsistencies between lifetime and polarizability measurements in Cs 
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Electric-dipole matrix elements for 6p — nd, n = 5, 6, 7 transitions in cesium are calculated using a 
relativistic all-order method. The resulting matrix elements are used to evaluate 5d lifetimes and 6p 
polarizabilities. The data are compared with experimental lifetime and polarizability measurements 
made by different groups. Domination of the 6p scalar polarizabilities by 5d — 6p dipole matrix 
elements facilitates an exacting consistency check of 5d lifetime and 6p polarizability data. Values 
of 5d — 6p matrix elements obtained from experimental 5d lifetime data are found to be inconsistent 
with those inferred from 6p polarizabilities derived from experimental Stark shift data. Our ab initio 
calculated 6p polarizabilities agree well with experimental determinations. 

PACS numbers: 31.15.Ar, 32.70.Cs, 32.10.Dk, 31.15.Dv 
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The understanding of the accuracy of ab initio calcu- 
lations in cesium is vital for the analysis of the Cs parity 
nonconservation (PNC) experiment I n 1999, moti- 
vated by a number of recent high- precision experiments, 
Bennett and Wieman reanalyzed the agreement of 
theoretical calculations and experimental data for a num- 
ber of Cs atomic properties and reduced the previous the- 
oretical uncertainty in the PNC amplitude by a factor of 
two. Utilizing measurements of the tensor transition po- 
larizability, (3, reported in same work, they demonstrated 
a 2.5a discrepancy between the value of the weak charge 
Qw predicted by the Standard Model and that derived 
from the Cs PNC experiment. Although several papers 
(for example, H S II 0, S @ ) , have addressed this 
disagreement since 1999, the issue of the accuracy of ab 
initio calculations in Cs continues to be of interest. 

In this work, we investigate the radiative properties 
of Cs 6p — nd transitions. Although these do not bear 
directly on PNC experiments done to date, they have 
been the subject of careful experimental investigation, 
and thus provide benchmarks for precise comparison of 
theory and experiment. In particular, there exist two in- 
dependent measurements of the lifetimes of the bd states 
fiol fTT| . which do not agree within their stated uncer- 
tainties. There also exist several experimental determi- 
nations of the 6p — 6s Stark shifts which allow to infer 
the values of polarizabilities of the 6p states^, HH. Il4j . 
Here we show that ab initio theory can check the mu- 
tual consistency of bd lifetime and 6p polarizability data, 
with an accuracy of about 1%. We find the lifetime and 
polarizability results to be inconsistent at this level. Our 
calculations agree with the experimental values of 6p po- 
larizabilities, but deviate from both determinations of 
the bd lifetimes. We suggest that further experiments 
are desirable in order to clarify this issue. In addition, 



understanding of the accuracy of the bd state properties 
in Cs is germane to the ongoing PNC experiment in iso- 
electronic Ba + , since the bd state is directly involved 
in this experiment. 

In outline, our approach uses a relativistic all-order 
method to calculate electric-dipole matrix elements for 
Cs 6p — nd transitions for n = 5, 6, and 7. These are 
used to evaluate bd radiative lifetimes and 6p polarizabil- 
ities (for the latter, we also include contributions from all 
other relevant states). Our calculations of the 6p scalar 
polarizabilities, which are in good agreement with exper- 
iment, show that they are dominated by contributions 
from bd — 6p transitions. These are the only electric- 
dipole transitions contributing to the bd state lifetimes 
(as we mention below, the bd — 6s electric quadrupole 
transition rates are negligibly small). Thus, it is possible 
to check consistency between polarizability and lifetime 
measurements by deriving bd — 6p matrix elements from 
bd lifetime measurements and substituting these values 
into the 6p polarizability calculations. For either of the 
two experimental lifetimes, 0, 0] this procedure yields 
a result that disagrees with directly measured polariz- 
abilities 

dill 

ll l| by several standard deviations. 
The particular all-order method used here is the lin- 
earized coupled-cluster method which sums infinite sets 
of many-body perturbation theory terms. We refer the 
reader to Refs. |l6l Il7l Il8| for a detailed description of 
the approach. The wave function of the valence electron 
v is represented as an expansion 
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where is the lowest-order atomic state function, which 
is taken to be the frozen-core Dirac-Hartree-Fock (DHF) 
wave function of a state v. This lowest-order atomic 
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TABLE I: Absolute values of electric-dipole 5d — 6p reduced matrix elements in Cs calculated in different approximations: 
Dirac-Hartree-Fock (DHF), third-order many-body perturbation theory (III), single-double all-order method (SD), single- 
double all-order method including partial triple contributions (SDpT) and the corresponding scaled values. R is the ratio of 
the 5d 3 /2 — to 5^3/2 — Gpi/2 transition matrix elements. All values are given in atomic units (edo, where ao is the Bohr 

radius) . 



Transition 


DHF 


III 


SD 


SD SC 


SDpT 


SDpT sc 


5d 3/2 - 6pi/ 2 


8.9784 


6.9231 


6.5809 


7.0634 


6.9103 


7.0127 


5d 3/2 - 6f> 3 /2 


4.0625 


3.1191 


2.9575 


3.1871 


3.1112 


3.1614 


R 


0.4525 


0.4505 


0.4494 


0.4512 


0.4502 


0.4508 


5d 5/2 - 6p 3/2 


12.1865 


9.4545 


9.0238 


9.6588 


9.4541 


9.5906 



state function can be written as \& v ) = al\0c), where 
0c) represent DHF wave function of a closed core. In 
equation a\ and are creation and annihilation op- 
erators, respectively. The indices m, n, and r designate 
excited states and indices a and b designate core states. 
The excitation coefficients p ma , p mv , p m nab, and p m nva 
are used to calculate matrix elements, which can be ex- 
pressed in the framework of the all-order method as lin- 
ear or quadratic functions of the excitation coefficients. 
We restrict the expansion given by Eq. to single and 
double (SD) excitations, with partial inclusion of triple 
excitations. The results obtained using the SD expan- 
sion are referred to as SD data throughout the paper 
and results obtained with partial addition of the triple 
excitations are referred to as SDpT data. We also per- 
formed third-order many-body perturbation theory cal- 
culations, following Ref. [l^, to better understand the 
size of higher-order correlation corrections. Unless stated 
otherwise, all results in this paper are expressed in the 
familiar system of atomic units, a.u., in which unit values 
are assigned to the elementary charge, e, the mass of the 
electron, m, and the reduced Planck constant h. 

Table Q] lists the 5c? — 6p reduced electric-dipole matrix 
elements in Cs as calculated using the Dirac-Hartree- 
Fock approximation (DHF), third-order perturbation 
theory (III), single-double all-order method (SD), and 
single-double all-order method including partial triple 
contributions (SDpT). We use semi-empirical scaling de- 
scribed, for example, in Ref. 01 to estimate some classes 
of the omitted high-order corrections. The scaled values 
are listed in rows labeled SD SC and SDpT sc . 

We use the 5c? — 6p matrix elements from Table [|] to 
calculate the lifetimes of the 5g? 3 / 2 and 5c? 5 / 2 levels in Cs. 
The Einstein A-coefficients A vw are calculated using the 
formula 0] 

_ 2.02613 xl0 15 |H|J>|M| 2 _ x 
vw ~ A3 2j v + l ' () 

where (u||.D||iu) is the reduced electric-dipole matrix el- 
ement for the transition between states v and w and A 
is corresponding wavelength in nm. The lifetime of the 
state v is calculated as 



The results are listed in Table [H] The experimental en- 
ergies from (2^ are used. The scaled SD values are taken 
as final values based on the comparison of a number of 
Rb, Cs, and Fr results jU 13] with experiment. The 
theoretical values differ substantially, by over 5%, from 
the experimental results (we note that that the exper- 
imental values from Refs. [13, El differ by 4%, which 
exceeds their stated uncertainties of 0.7% and 1%, re- 
spectively). One possible source of such a discrepancy is 
the contribution of the 5d — 6s electric-quadrupole tran- 
sition to the 5c? lifetime. Our calculation of this rate, 
using the all-order method, yields a corresponding Ein- 
stein A-coefficient for the 5g? 5 / 2 — 6s transition of 19 Hz, 
which is only 0.02% of the corresponding electric-dipole 
A-coefficient of 741 kHz (see Table [0J. Thus, the contri- 
bution of the electric-quadrupole transition to 5c? lifetime 
is entirely negligible within the present experimental and 
theoretical uncertainties. 

To clarify such a large disagreement we check the con- 
sistency of the experimental 5c? lifetime measurements 
with Qp polarizability measurements, which involves con- 
tributions from the same transitions. First, we use ex- 
perimental 5c? lifetimes from 10] to determine the 5c?— Qp 
reduced matrix elements. Inverting Eq. J^J), we find for 
the 5g?5/ 2 — 6p 3 /2 matrix element: 

\(5d 5/2 \\D\\6p 3/2 }\= 9.916(35). (4) 

To derive the 5c? 3/ / 2 — &Pi/2 an d 5c? 3 / 2 — 6p 3 / 2 matrix ele- 
ment, the lifetime of the 5c? 3 / 2 level alone is not sufficient 
and some assumption about the ratio R of these matrix 
elements must be made. We use the theoretical SD SC 
value 0.4512(18) from Table [I] for the ratio and assume 
the deviation of other high-precision theoretical results 
in Table [I] from this value to be its uncertainty. The vari- 
ation of the ratio from one approximation to another is 
far smaller than the variation in the individual matrix 
elements, thus the uncertainty is rather low (0.4%). The 
resulting values of the 5c? 3 / 2 — 6p matrix elements are: 

|<5c? 3/2 p||6p 1/2 )| = 7.283(60), 
K5c? 3/2 p||6p 3/2 )| = 3.286(27)(13). (5) 

We separated the uncertainties in the 5c? 3 / 2 — 6p 3 / 2 matrix 
elements into contributions from the 5c? 3 / 2 lifetime mea- 
surement (0.027) and from the estimation of R (0.013). 
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TABLE II: The values of Einstein A-cocfficients A vw (in MHz) and final lifetimes (in ns) for 5ds/2 and 5^3/2 states in Cs. The 
theoretical values are compared with experimental results from [TfH ] and [Til ]. 



Level 


Transition 




SD 


SD SC 


SDpT 


SDpT sc 


Expt. [101 


Expt. [Ill 


5d 5 /2 


5d 5/2 - 6p 3/ 2 


A-vw 


0.646 


0.741 


0.710 


0.730 










r 


1547 


1350 


1409 


1369 


1281(9) 


1226(12) 


5d 3/2 


5d 3/2 - 6pi /2 




0.804 


0.926 


0.886 


0.913 






5d 3/2 - 6p 3/2 


^™ 


0.094 


0.109 


0.104 


0.107 








T 


1114 


966 


1010 


981 


909(15) 





TABLE III: Contributions to the 6pi/2 and 6P3/2 scalar po- 
larizabilities ao in Cs and their uncertainties Sao, in units of 
ag. The values of corresponding matrix elements d (in a.u.), 
their sources and uncertainties Sd (in %) are also given. The 
6p — 6d and 6p — Id matrix elements are from the present 
SDpT all-order calculation. 



Q (6p 


1/2) 


d 


Sd 




QO 


Sao 


6f>l/2 - 


5^3/2 


-7.283 


0.8 


l-i-U] 


1168.4 


18.7 


6P1/2 - 


6s 


-4.489 


0.1 


[20] 


-131.9 


-0.3 


6f>l/2 - 


6(^3/2 


4.145 


4.8 


SDpT 


110.2 


10.6 


6P1/2 - 


7s 


-4.236 


0.5 


[6] 


178.4 


1.8 


6pi/2 - 


7d 3 / 2 


2.033 


1.7 


SDpT 


20.3 


0.7 


6P1/2 - 


8s 


-1.026 


0.6 


[6] 


5.9 


0.1 


6P1/2 - 


9.s 


0.550 


0.5 


m 


1.4 


0.0 


«tail 








DHF 


35.4 


10.6 


Q^corc 








[21J 


15.8 


0.3 


Total 










1404 


24 


ao(6p3/ 2 ) 


d 


Sd 




ao 


8a 


6P3/2 - 


5^3/2 


3.286 


0.9 


no] a 


142.7 


2.6 


6P3/2 - 


5d 5 /2 


9.916 


0.3 


[10] 


1255.5 


8.8 


6P3/2 - 


6s 


-6.324 


0.1 


[20] 


-124.7 


-0.2 


6P3/2 - 


6^3/2 


-2.053 


4.6 


SDpT 


14.2 


1.3 


6P3/2 - 


6^5/2 


-6.010 


4.3 


SDpT 


121.2 


10.4 


6P3/2 - 


7s 


-6.473 


0.5 


[6j 


225.3 


2.3 


6P3/2 - 


7d 3 / 2 


-0.969 


1.5 


SDpT 


2.4 


0.1 


6P3/2 - 


7d 5/2 


-2.868 


1.4 


SDpT 


21.0 


0.6 


6P3/2 - 


8s 


-1.462 


0.6 




6.2 


0.1 


6P3/2 - 


9s 


0.774 


0.6 


m 


1.4 


0.0 


Qtail 








DHF 


38.7 


11.6 


Q^corc 










15.8 


0.3 


Total 










1720 


18 



a Derived from the experimental 5d 3 j2 lifetime ll(l using theoret- 
ical ratio of the 6p 3 /2 — ^3/2 an< l §Pl/2 ~ ^3/2 matrix elements. 



Combining them, we obtain 3.286(30). The contribution 
of the uncertainty in R to the uncertainty in the value of 
5^3/2 — 6j?i/2 matrix element is negligible. 

The scalar ao and tensor 012 polarizabilities of of an 
atomic state v are calculated using formulas 



, , E n — E v 



n \\ D \W 



= 4 



(6) 



6(j v + l){2j v + 1)(2 A 

jv+jn + 1 ) 3v 1 jn 

1 jv 2 



1/2 



3), 
(n\\D\\v) 2 

E U — Ey 



(7) 



TABLE IV: Contributions to the 6p 3 / 2 tensor polarizability 
«2 in Cs and their uncertainties Sa2 in a^. The values of 
corresponding matrix elements d (in a.u.), their sources and 
uncertainties Sd (in %) also given. 







d 


Sd 




Q2 


5a2 


6P3/2 


- 5d 3 / 2 


3.286 


0.9 


[10j a 


114.2 


2.1 


6P3/2 


- 5<f 5/2 


9.916 


0.3 


no] 


-251.1 


-1.8 


6P3/2 


-6s 


-6.324 


0.1 


[20] 


124.7 


0.2 


6P3/2 


- 6d 3/2 


-2.053 


4.6 


SDpT 


11.4 


1.0 


6P3/2 


- 6d 5 / 2 


-6.010 


4.3 


SDpT 


-24.2 


-2.1 


6P3/2 


- 7s 


-6.473 


0.5 


m 


-225.3 


-2.3 


6P3/2 


- 7d 3/2 


-0.969 


1.5 


SDpT 


1.9 


0.1 


6P3/2 


- 7d 5/2 


-2.868 


1.4 


SDpT 


-4.2 


-0.1 


6P3/2 


-8s 


-1.462 


0.6 


M 


-6.2 


-0.1 


6P3/2 


-9s 


0.774 


0.6 


M 


-1.4 


0.0 


Oxtail 








DHF 


-7.0 


2.1 


Total 










-267.3 


4.7 



"Derived from the experimental 5<i 3 /2 lifetime using theoret- 
ical ratio of the 6P3/2 — 5^3/2 an< l &Pi/2 ~ 5^3/2 matrix elements. 



where D is the dipole operator and formula for ao in- 
cludes only valence part of the polarizability. The main 
contributions to the polarizability, a ma j n , come from 
transitions between 6s, 7s, 8s, 9s, 6p, 7p, Sp, 9p, 5d, 
6d, and 7d levels; the remainder, ataii, is calculated 
from summing over all other valence-excited states of 
the system (which is confined in a sphere of radius 
75 ao). The core contribution to the scalar polarizability 



15.8 0%, 



is taken from 21], where it was cal- 



culated in random- phase-approximation (RPA). We note 
that this value includes the contribution from the valence 
shell and, therefore, must be compensated by the addi- 
tional term, a vc , which is equal to the contribution from 
the valence shell divided by (2j v + 1) with an opposite 
sign. We find that the a vc term is negligible for np states 
and very small (below 0.2%) for the 6s state. We list the 
contributions to Cs 6p polarizabilities in Tables IIIII and 
IIVI The corresponding electric-dipole matrix elements d, 
their sources, and uncertainties Sd are also given. The 
values for (is — np and 7s — np transitions are taken from 
Ref. [6| , where the "best value" set of these matrix ele- 
ments was compiled for the calculation of the tensor tran- 
sition polarizability (3. The 6p — 6d and 6p — 7d matrix 
elements are from the present ab initio SDpT calculation. 
The values of the 5d — 6p matrix elements are derived 
from the 5d lifetime experiment [lOj. The same data set 
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TABLE V: Calculated and experimental values of Cs polarizabilities, in aj). Calculation (a) uses 5d — 6p matrix elements data 
derived from the 5d lifetime experiment ^jj (results of Tables ITTT1 and IT\^ : calculation (b) uses 5d — 6p theoretical all-order 
values (SD scaled data). All other contributions in calculations (a) and (b) are the same. 







(a) 

Expt. 5d — 6p 


Present 

(b) 

Theory 5d — 6p 


Expt. [12] 


Expt. [13] 


Expt. [14] 


Qo(6p3/ 2 ) 


— ao(6s) 


1322(18) 


1248 


1264(13) 


1240.2(24) 




«o(6pi /2 ) 


- ao(6s) 


1006(24) 


936 


970(9) 




927.35(12) 


a2(6p 3/2 ) 




-267(4.7) 


-261.2 


-261(8) 


-262.4(15) 





is used in both tables. The uncertainties of all contri- 
butions are listed separately. The uncertainties listed in 
Ref. [|| are used for 6s — np and 7s — np transitions; the 
difference between SD and SDpT data is taken to be the 
uncertainty of the 6p — 6c? and 6p — 7d matrix elements 
calculated in this work. The uncertainties of the 5d — 6p 
matrix elements obtained from the lifetime experiment 
|10| are derived above. The uncertainty of the core term 
dcore is taken to be 2% based on the comparison of RPA 
data for closed core systems with experiments and high- 
precision calculations. The uncertainty of the remaining 
contribution a ta ii is estimated to be 30% based on the 
comparison of the DHF results with correlated values. 

We also calculate the scalar polarizability of the 6s 
state using the same methods and data set as for the 6p 
polarizability. The resulting value ao(6s) = 398.2(0.9) Oq 
and its uncertainty are dominated by contributions of 
the 6s — 6p matrix element taken from experiment of 
Ref. [2(j ■ We use this result when calculating differences 
of 6p and 6s polarizabilities. The recent measurement 
of the ground state polarizability in Cs yielded the value 
a (6s) = 401.0(0.6) ag [25|. 

We compare the final results for the differences of 
the 6p and 6s scalar polarizabilities ao and the ten- 
sor polarizability a 2 with experiment in Table [3 The 
results of the above calculation (data from Table IIIII 
I1V|) . where we used 5c? — 6p matrix elements derived 
from the 5c? lifetime experiment are listed in column 
(a). We find that the difference of the 6p3/ 2 and 6s 
scalar polarizabilities which uses numbers for 5c? — 6p 
matrix elements derived from |10| 5c? lifetime measure- 
ments ao(6p 3 / 2 ) — cto(6s) = 1322(18) afj is inconsistent 
with both experimental values 1240.2(24) ag and 
1264(13) al The difference with first value, which 

has the smaller uncertainty is 4.5cr and the difference with 
the second value is 2.6c. The difference of the &P1/2 and 
6s scalar polarizabilities which uses numbers for 5c? — 6p 
matrix elements derived from ^(j 5c? lifetime measure- 
ments ao(6pi/ 2 ) — ao(6s) = 1006(24) ag is also inconsis- 
tent with the most recent and most precise experimental 



value 927.35(12) a% pj by 3.2cr. The value for the 6p 3/2 
tensor polarizability —267.3(4.7) Oq has much larger un- 
certainty owing to strong cancellation of the contribu- 
tions from different transitions, and the difference is la. 
We note that if we were to use another 5g? 5 / 2 lifetime ex- 
periment [TTj . the discrepancies with polarizability mea- 
surements onl y in crease. Thus, neither 5dn/? lifetime ex- 
periment 0, ^3 is consistent with either |12| , p| , or 
|14| Stark shift measurements within the quoted uncer- 
tainties. 

We calculate that the experimental value of 
ao(6pi/ 2 ) - a (6s) = 927.35(12) ag 1 4] corresponds to 
the lifetime of the 5g?3/ 2 state Tsd 3/2 = 975(14) ns and the 
experimental value of ao(6p3/ 2 ) — ao (6s) = 1240.2(2.4) 
|13| corresponds to the lifetime of the 5c? 5 / 2 state 
T 5d 5/2 = 1359(18) ns. The uncertainties in these life- 
time values are dominated by the uncertainties in the 
values of 6p — 6c? transitions and the uncertainty in the 
contribution cttaii as evident from Table [TTT1 

Finally, we repeated the polarizability calculation by 
replacing the 5c? — 6p matrix elements derived from the 
lifetime experiment by our theoretical values (SD SC ) from 
Table |U All other matrix elements and contributions are 
exactly the same as in the first calculation. The results 
are listed in column (b) of Table As expected, they 
are quite different from the previous calculation (a) as 
our theoretical 5c? — 6p matrix elements are substantially 
different from the values derived from 5c? lifetimes. We 
find that our theoretical polarizability data are in good 
agreement (0.4%-l%) with experimental results. 

In conclusion, we find the experimental measurements 
of 5c? lifetime and 6p scalar polarizabilities to be inconsis- 
tent within the uncertainties quoted by the experimental 
groups. Our theoretical calculations are consistent with 
polarizability experiments but not with the lifetime mea- 
surements. Thus, further measurements of the properties 
of 5g! and 6p states are of great interest for clarification of 
this issue and for providing benchmark values for 5o? — Qp 
matrix elements. 
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